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Abstract: We consider a simple extension of the Standard Model by the addition of N 
real scalar gauge singlets (p that are candidates for Dark Matter. By collecting theoretical 
and experimental constraints we determine the space of allowed parameters of the model. 
The possibility of ameliorating the little hierarchy problem within the multisinglet model 
is discussed. The Spergel-Steinhardt solution of the Dark Matter density cusp problem 
is revisited. It is shown that fitting the recent CRESST-II data for Dark Matter nucleus 
scattering implies that the standard Higgs boson decays predominantly into pairs of Dark 
Matter scalars. It that case discovery of the Higgs boson at LHC and Tevatron is impossible. 
The most hkely mass of the dark scalars is in the range 15 GeV < < 50 GeV with 
BR{h (pip) up to 96%. 
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1 Introduction 

The evidence for dark matter (DM) has now become almost overwhelming [1], yet the 
standard model (SM) does not contain a viable DM candidate, which necessitates a modi- 
fication of this model. The simplest possible such extension consists of introducing a real 
scalar field which is a gauge singlet. Many authors considered this singlet extension of the 
SM; the earliest publication we are aware of was by Veltman and Yndurain in [2] (though 
their motivation was different). The DM issue was first addressed in this context in [3] and 
[4] and then followed by other authors [5]-[15]. 

In this paper we will consider a generalization of this model by extending the DM sector 
to an unbroken (global) 0(N) model. We will derive relevant theoretical and experimental 
constraints on this model and determine range of allowed parameters. Since we demand the 
new singlets do not acquire a vacuum expectation value, the effects of the dark sector on 
precision observables is much suppressed, appearing only at the two (or higher) loop level. 
The only experimental constraints are then obtained from the requirement that the scalars 
are adequate candidates for DM. The theoretical constraints on the model are derived from 
vacuum stability, perturbativity, and triviality. In addition we will also determine the fine 
tuning conditions that ameliorate the little hierarchy problem by seeking model parameters 
that tame quadratically divergent contributions to the Higgs boson mass [16], [17] following 
strategies adopted in [18]. 
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The model discussed below is also interesting in light of the DM signals recently an- 
nounced by the CRESST-II collaboration [19]. We will see that, not only the most likely 
regions found by CRESST-II are consistent with the model, but in a region of parameter 
space where the SM Higgs decays invisibly into DM pairs. Such a Higgs boson would not 
be detected at either the Tevatron or LHC. 

The paper is organized as follows. In Sec. 2 we introduce the model and formulate 
strategy adopted in this paper for its study. In Sec. 3 we review the vacuum stability 
and unitarity bounds that restrict the parameter space of the model, while in Sec. 4 the 
issue of triviality is addressed. Sec. 5 is devoted to a discussion of one- and leading two- 
loop corrections to the Higgs boson mass in the presence of extra singlets. In Sec. 6 the 
present Dark Matter abundance is investigated for the Cold Dark Matter (CDM) and 
Feebly Interacting Dark Matter (FIDM) scenarios. Sec. 7 shows results of constraining 
the model by results of direct searches of DM. In Sec. 8 we revisit the Spergel-Steinhardt 
hypothesis of self-interacting DM. Sec. 9 contains summary and conclusions. 

2 The Model 

We consider the SM of electroweak interactions extended by the addition of N scalars (f that 
are singlets under the SM gauge group SU{3) x SU{2) x ?7(1), and transforming according 
to the fundamental representation of 0{N), under which all SM fields are singlets; for 
simplicity we assume that 0{N) is an exact symmetry of the model ^. With the intention 
of providing a Dark Matter (DM) candidate we impose an additional Z2 symmetry, which 
is not spontaneously broken, and under which (f is odd: (f —0 while all other fields are 
even ^ . The most general, symmetric and renormalizable potential reads: 

V{H, ^) = -^i^H^H + XniH^Hf + ^/.^ ^ (^2^2 ^ ^^^t^^2 ^ ^2.1) 

where H is the SM SU{2) Higgs isodoublet. The Lagrangian density for the scalar sector 
is then given by: 

Lscaiar = + D^H^D^'H - V{H, 0) . (2.2) 

As usually the minimum of the potential breaks spontaneously electroweak symmetry via 
non-zero vacuum expectation value of the Higgs doublet (H) = (0,'u/\/2), v = 246 GeV. 
Since we require the 0{N) symmetry {Z2 for = 1) to remain unbroken, we assume 
that fj!^ > 0, so {(f) = 0. Note that {(f) = implies no mass-mixing between (f and 
i7, therefore the existing collider limits on the Higgs properties are not modified. After 
the symmetry breaking the physical scalars have masses = "^1/ + 3Xhv^ — 2jj!jj and 
= + ^x'^^- Note that all components of (f have the same mass as the consequence 
of 0{N), this degeneracy can be removed by adding a generic mass term (fi^)ij(fi(fj that 
only breaks 0{N) softly. 

^This condition could be relaxed without significantly altering the conclusions; the only consequence 
would be the presence of many more parameters of the model. 

^For N > 1 this Z2 symmetry is a consequence of 0{N) invariance. 
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The model then contains five unknown parameters: rrih^m^, Xx, X(^, N . Our goal is 
to constrain the parameters taking into account available restrictions: theoretical (vacuum 
stability, unitarity/perturbativity, triviality of the scalar sector, Higgs mass correction fine- 
tuning) and experimental (DM relic abundance, direct detection experiments). 

3 Stability and perturbative unitarity 

In order to stabilize vacuum we will assume that the scalar potential (2.1) is (at tree level) 
bounded from below. This condition implies 



A(^Aif _ rrih I X^ 



Xh,X,>0; A.>-y^ = -^yf. (3.1) 

Tree-level unitarity constraints emerge from the SM condition for VlVl scattering [20] 
and from the requirement that all possible scalar-scalar scattering amplitudes are consistent 
with unitarity of the S matrix [21] 

ml < — v'^, X^ < Stt and \Xx\ < 47r . (3.2) 
3 

Finally, the condition that the global 0{N) symmetry remains unbroken requires /x^ > 
which leads to the very useful inequality: 

ml > Xxv^ ; (3.3) 

a consequence is that light scalars (m^ <C v) must couple very weakly to the SM (A^; <^ 1). 

4 Triviality 

The 'triviality bound' is a constraint on the SM generated by the requirement that under 
the renormalization group evolution the quartic coupling constant Xh remains finite up to 
the UV cut-off scale A of the model; in other words, that the Higgs-boson Landau pole 
should be above A. This requirement implies conditions on the initial values of various 
running parameters of the model, and in particular on Xh{i^ = ^W"); this, in turn, leads 
to an upper limit on the Higgs boson mass as a function of A. 

In order to determine the location of the Landau pole as a function of the Higgs mass 
one has to solve the renormalization group evolution (RGE) equations for all of the running 
parameters of the model. We will be interested only in cutoff scales below 50 TeV, in which 
case the RGE of the gauge and top-quark Yukawa couplings can be safely neglected. Thus 
we only need to consider the evolution of Xh, X^ and A^; as determined by 

levrV^ = Igf + Igl + \9hl - ^vt + 24a1, + 12^^^ - 35? - '^glXn 

+2N\l (4.1) 

2 dXx ^ A .^9 39 9 



IGttV^ = A, \12\h + A^ + 8A, + 6y/ - -gi - -g^j (4.2) 
167rV^ = 48A2 + l(8 + iV)A2 (4.3) 
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where ^i, ^3 are the gauge couphng for U{1)^ SU(2)^ SU{3)^ respectively, and yt is the 
top quark Yukawa couphng Unique solutions could be obtained once the initial conditions 



(4.4) 
(4.5) 
(4.6) 



are specified. 

For a given A we shall require that there is no pole in the evolution of scalar quartic 
coupling constants at energies below A. The region in the (A, rrih — v\/2Xho) plane allowed 
by this constraint, depends on the initial parameters A^^o, ^^pO, and the number of scalars 
N. For each A the maximum allowed value of ruh constitutes the triviality bound on the 
Higgs mass. One might worry that each quartic coupling constant will have a different pole 
location, this, however, does not occur when A^^o 7^ (as we assume), for then all the RGE 
equations are coupled and \h,lp,x diverge at the same value of ii. 

The triviality bound for mh as a function of A is illustrated in Fig. 1. Note that the 
allowed region shrinks as A^^o grows for fixed N ^ and as N increases for fixed A^^o, lowering 
the upper bound on ruh in either case. This behavior is a direct consequence of the term 
2N\^ in (4.1): increasing N and/or A^; amplifies the evolution of \h- For a given (A,m/i) 
and fixed A^^o there is a range of A^^o foi" which the Landau pole occurs above A, as shown 
in Fig. 3. Note the asymmetry of the allowed (inner) region, which is a consequence of the 
8A2 term in (4.2). 

In Fig. 2 we show triviality limits for A^^ as a function of A for the case A^^o = when 
the evolution decouples from the SM. 

5 Fine-tuning in Higgs mass corrections 

If we assume that the SM is an effective theory valid at scales below the UV cutoff A, 
radiative corrections [22] shift the Higgs mass to a scale of order A/\/47r. The existing 
bounds on the scale of new physics imply A > 0( TeV), which is in conflict with a Higgs 
mass below 1 TeV range as required by the unitarity constraint [20] . A way of alleviating 
this conflict is to arrange cancellations to occur within the 0{h?) contributions to m^. 
These cancellations can be obtained through a fine tuning of the parameters or through 
a an extension of the SM together with an appropriate new symmetry (as is the case in 
supersymmetric models). 

The first solution to this problem was proposed by Veltman [22] who observed that 
the 1-loop quadratic corrections to rrih would vanish if + m\ + 2m^ — 4m^ = 0, but 
this possibility is now experimentally disallowed. Here we shall follow a similar approach 
within the multi-singlet extension of the SM. 

In pursuing this approach it is important to determine the possible effects generated 
by higher order radiative corrections. In a theory with many couplings A^ the general form 
of leading higher-order contributions (those containing the highest power of log (A)) to the 
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Figure 1. The triviality upper bound on the Higgs mass as a function of the cut-off A for X^{mw) = 
0.1, for A^^o = 0.1,1,2 (black, blue and red curves, respectively) and N = 1,6,12 (starting with 
the uppermost curve); the region above each curve is excluded the by the triviality constraint for 
the corresponding set of parameters. For A^^o = 0-1 non-standard effects are so small that upper 
boundaries for = 1, 6 and 12 are indistinguishable from the SM, accordingly the black curves are 
nearly identical to the bound derived within the SM. 
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Figure 2. Maximum A^^o allowed by the triviality condition as a function of A for N = 1, 6, 12 and 
^xo =0- The right (left) panel is for initial conditions provided at a 10 MeV (100 MeV) 



quadratically divergent contributions to at n + 1 loops take the form [23] 



n=0 



log 



(5.1) 
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Figure 3. Maximum and minimum X^o allowed by the triviality bound for A = 10 TeV, A^g = 0.1 
and N = 1,6,12. 



where fi is the renormahzation scale and the coefficients satisfy 

d_ 

'dXi- 



{n + l)fn+i 



^ f 



(5.2) 



Thus in order to determine the two-loop leading logarithmic corrections to the Higgs 
mass one needs the one-loop corrections and the first order beta functions (for generic 
results and applications to the single scalar case see [24]). The results for A^-singlet model 
read: 

A2 



^''^hl-loop 



167r2 V 

1 



12yf + ^gf + \9l 



5m 



167r2 
A2 



h 2— loops 



2\2 



6XHml log 

log (-) [25gt + 9glgl - 15gl + 34^?^^ ^ 54^2^2 



2 

mt 



+ 2X^711 log 



+ A^ 



(5.3) 



(167r 

+192^2^2 _ imytXH - 365? - + lUy^XH + 288A1^ 

-^NglX^ - mglX^ + 12Ny'^X^ + 24iVAijA^ + 40iVA2 + 2iVA^A<^](5.4) 

the logarithmic terms (oc m^ log A) in the one- loop correction were kept since they are 
relevant in the range mh <C m^ < A) ; terms oc log A are always numerically negligible 
and are included for completeness. For the numerical results the renormahzation scale was 
chosen to be the vacuum expectation value of the Higgs field, /i — 246 GeV. The SM result 



can be recovered in the limit A^ 



A. 



0. 



We implement the fine tuning condition in a standard manner [18] by requiring 

< A;, (5.5) 



5ml 




^Kl-loop + ^^L-loops 
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and determine the region in the (A, rrih) plane where this is obeyed for a given choice of 
the fine tuning parameter Ah. Specificahy, for each choice of (A,m/j) with fixed N and 
Xx we check if there exist and A^^ such that condition (5.5) is satisfied; the results are 
presented in Figs. 4 and 5. 

It is worth noting that the two-loop (f corrections to ruh (5.4) grow as Xx/ (IGtt^) ln(A//x) 
relative to the one-loop contributions(5.3). Perturbative consistency then requires this fac- 
tor to be < 1, which translates into a maximal cutoff A^ax beyond which perturbation 
theory is no valid: 

A < An.ax = (5.6) 

Therefore the regions of large A in Figs. 4 and 5 should be considered with a certain caution 
(that should increase with A^;). The range of A is extended up to 10^ GeV only for the 
purpose of the SM (A^; = 0) and small A^; cases (A^; < 0.95 for = 246 GeV). 

We now consider the regions in the (A, ruh) excluded by the triviality and fine tuning 
constraints, first for A^^ > and then for A^^ < 0. 

Xx ^ 0. As one can see in Fig. 4, the fine-tuning condition (5.5) defines a region in the 
(A, rrih) plane with a characteristic funnel-shaped boundary (bounded, for various values 
of Xx, by the solid lines shown in the figure). Given the experimental bounds [25] on m/^, 
115 GeV < rrih < 141 GeV, we see that as A^; or N grows, the lower branch of the funnel 
allows larger regions of A for fixed A, or smaller values of Ah for fixed A; in contrast the 
upper branch of the funnel allows smaller regions of A for fixed A/^, or larger values of Ah 
for fixed A. As a result the combined effect is complicated, for example the region allowed 
by (5.5) and the experimental limits on ruh shrinks as A^; goes from 1.25 to 2 when = 1, 
and as A^; goes from 0.5 to 1 for = 3. For fixed A^^ and N the allowed range of A shrinks 
as Ah becomes smaller. 

Several comments are in order here. For the SM (A^^ = 0) there is a disallowed region 
for A between 1-10 TeV and - 10^ TeV; this region shrinks as A grows. For instance, if 
A 10^ GeV then the present limits on the Higgs boson mass, can be accommodated only 
for large fine tuning: Ah > 100. In order to have Ah < 100 a modification of the SM must 
be introduced, or low (^ 1 TeV) or very high > 10^ TeV values of A must be assumed. 

Note that for each point (A, ruh) below the SM (A^; = 0) upper branch there exist A^; 
such, that the point is allowed. Therefore in the model we discuss here, for a fixed cutoff 
A, much smaller Higgs-boson masses are allowed. In particular the range of 115 GeV < 
ruh < 141 GeV could be easily accommodated with a cutoff much larger than the SM one. 

Let's concentrate on the cutoff A = 10^ GeV. As it is seen from the left panels is 
Fig. 4 for such A, in order to have the range of presently allowed Higgs masses 115 GeV < 
rrih < 141 GeV within the region of fine tuning 6mf^/mf^ below 10, one needs A^; of the 
order of 1.2 — 1.3 and 0.5 — 0.6 for = 1 and = 3, respectively. However, if we relax 
the fine tuning condition so that 5m\/m\ < 100 (see the right panels of Fig. 4) then, as 
expected, lower A^; is sufficient; Xx — I for 115 GeV < ruh < 130 GeV and A^; 2^ 0.5 for 
130 GeV <mh< 141 GeV, for A^ = 1. If A^ = 3, then A^: ^ 0.5 is large enough to bring the 
whole ruh range into the region of the fine tuning 6mf^/mf^ < 100. In general, A^; = 0{1) 
is needed to ameliorate the SM fine tuning. 
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Figure 4. Constraints on and A for Ah = 1, 10 and N = 1^3. The horizontal band shows the 
presently allowed range for the Higgs-boson mass [25]: 115 GeV < rrih < 141 GeV. Dashed curves 
correspond to the triviality limits for the SM (A^^ = 0, dashed black curve) and singlet extensions 
for several values of A^^ (dashed blue and red curves); regions below the curves are allowed by this 
constraint. The solid lines correspond to the two-loop fine-tuning constraint for the SM (A^^ = 0, 
solid black) and singlet extensions for several values of A^; (solid red and blue curves), the regions 
to the left of the corresponding funnels are allowed by this constraint. Every point inside the dark 
(light) gray region is forbidden by the triviality (find tuning) condition for any value of A^^. For every 
allowed point outside the dark gray and gray regions there exist |Aa^| < Stt, 10 GeV < m^^ < 10^ GeV 
and < A(^ < 2 such that (5.5) is fulfilled. 

On the other hand, one should remember that for each choice of m^, Xx is restricted 
from above by the condition (3.3), therefore if 0.5 < A^; < 1.5, then the corresponding (f 
mass must be in the range 174 GeV < < 301 GeV. Alternatively, if we require = 
15 — 60 GeV (the range relevant for direct DM detection, see Sec. 7) then 0.06 < A^; < 0.24; 
as seen from Fig. 4 such low (f masses are inconsistent with = 100, A = 10^ GeV and 
= 1 or 3. In order to allow for that low mass, either larger A^ > 100 or lower cutoff 
A would be necessary. 
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Figure 5. Similar as Fig. 4 for negative A; 



It is worth noting that for large ruh {rrih ^ 310 GeV) the SM contribution to the 
quadratic divergence is dominated by the Higgs boson, therefore in order to enlarge the 
allowed region the extra contribution from (f in the loop should be of the opposite sign 
(corresponding to A^; < 0). 

Xx < 0. In this case the limit (3.3) provides no restriction, on the other hand the vacuum 
stability condition (3.1) becomes important. This constraint allows |Aa;| > 1 but only for 
ruh > 170 GeV due to the unitarity limit (3.2); saturating this limit by taking A^^ = Stt 
requires |Aa^| < 0.65 — 0.85 for 110 GeV < rrih < 141 GeV. In this range the allowed values 
of A are more restricted when the singlets are included than for the SM, this is because 
for Xx < the 1-loop (^-contribution to 6mj^ is of the same sign as the top-quark one, 
therefore for small A^^ funnels are shifted towards larger values of m/^, as needed by (5.5) 
in order to compensate the top and (f effects; as A^; increases the funnels move downwards 
as a consequence of the growing 2-loop contributions to ^"^^2— loops' 

The dashed curves in Figs. 4 and 5 give the triviality limits on rrih for various values 
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of \x (see also Fig. 1). As seen from the figures, this constraint is not important for 
115 GeV <mh< 141 GeV. 

Summarizing this section, one can say, that for < A^; < 1 the model allows to shift 
the UV cutoff up to A < 10^ GeV keeping the Higgs-boson mass within the experimentally 
allowed range 115 GeV < ruh < 141 GeV. The case A^^ < is disfavored, as the cutoff 
in this case appears to be lower than in the SM as long as we stay within the range of 
parameters allowed by perturbative expansion. 



6 DM relic abundance 

The singlet being odd under the Z2 symmetry, is stable, and therefore a dark matter 
candidate, in which case the parameters of the theory should also lead to the DM abundance 
consistent with the WMAP observations, d^"^^ = 0.110±0.018 (3cr) [26]. In order to gauge 
the consequences of this constraint on our model we first calculate the relic DM abundance 
by solving the appropriate Boltzmann equation. Considering first the single component 
(TV = 1) case, we have [27]: 



dT^^^^ ^^^^^^^B^^ M^^^ J (2^)3exp(£;/T)±lf^) 

where / = n/T^ ^ n is the particle number density of the DM candidate (in our case the 
singlet), fEqiT) is the equilibrium distribution, ^^(T) is the number of relativistic degrees 
of freedom at temperature T, mpi is the Planck mass and {av) is the thermally-averaged 
cross section for DM + DM SM + SM annihilation processes ([28], [10]): 



i(s) = 2^s(s-4my,7(s) (6.3) 

where x — vn^lT ^ Ki^2{x) are the modified Bessel functions, and a{s) is the DM -\- DM 
SM + SM annihilation cross section normalized in the standard manner. In analytical 
estimates is useful to note that {av) ^ a(4m^)/4m^ for non-relativistic DM; however we 
do not use this approximation in the numerical results shown below. 

In the case of A^-singlet scalar DM candidates with 0{N) symmetry, all components 
contribute equally to the relic abundance, just like a single field with N degrees of freedom. 
Therefore the total abundance equals 

<M = E^DM = ^^f^DM (6.4) 

i 

where f^j^^ dark matter relic density from the z-th scalar field. By virtue of the 0{N) 

symmetry it is sufficient to consider the Boltzmann equation for one of the components of 

In the following we will discuss two different limiting cases for the Boltzmann equation: 
Cold Dark Matter (CDM) and Feebly Interacting Dark Matter (FIDM). We do not consider 
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Figure 6. Feynman diagrams illustrating (p(p annihilation into SM particles. 

the hot dark matter solution, since it is inconsistent with structure formation at the galaxy 
scale, see e.g. [29]. 

The diagrams contributing to 00 annihilation into SM particles are shown in Fig. 6. 
The corresponding cross sections (for = 1) are available in the literature (e.g. [5] and 
[10]); we have verified the results of [10]: 

aww{s) - ^ V 1 - , (s-ml)^ + mlVl ^ + ^ 



s {s - mlY + m^r^ 
<Jzz{s) = 7^\/l ^7 ^ 2t^2 (^2^ ^ + 1 




3 

4mj \ mj s 



where F(^) = ArcTanh(^)/^, ^ = y (5 - 4m^)(5 - 4m^)/(5-2m^). The total cross section 
is then 

a{s) = aww{s) + o-zz(<s) + ^ ^//l"^) + (6-6) 

/ 

where the sum runs over all fermions /. 

When > 1 one has to modify the Higgs width as it can decay into N possible 
final states: 

T/z ^h^SM + ^h^cpcp (6.7) 

Nv^X^ r-o 

^h^^^ = ^^^2 y^h - ^H{mh - 2m^) (6.8) 

where Oh is the Heaviside step function. 
6.1 Cold Dark Matter 

In the CDM approximation the solution to the Boltzmann equation for the 0{N) model 
yields the relic density [27] 

^N^h^ = = 1.06 X io9 J^'^f \ (6.9) 



Pcrit ' \/(hmpi (av) GeV 
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where Xf = m^/Tf^ and Tf is the freeze-out temperature given in the first approximation 

by 

Xf = log ( O.OSsi^^l^^^ (6.10) 

where {av) is the thermally averaged cross section for + SM + SM. Note that 

the A^-dependence of enters through the Higgs width (6.7) contained in {av)] this weak 
dependence is furthered softened by the log function in (6.10). Therefore it is reasonable 
to neglect the N dependence in (at least in the parameter range where the perturbative 
expansion is valid). Note also that in the case where the t- and ix-channel annihilation 
diagrams with (p exchange can be neglected (e.g. when m^p » ruh and/or if A^^ <^ 1), the 
abundance depends only on the combination A^A^. 

For a given choice of N and (rrih^m^) we will look for A^; such that the constraint 
f^J^^P) = 0.110 ± 0.018 is satisfied: 

0.092 < S^^M < 0-128. (6.11) 

Even though (6.9) often leads to accurate results for the abundance, in deriving our numer- 
ical results we used micrOMEGAs [30] - a code dedicated for calculation of DM properties. 
Our results are presented in Fig. 7, where we have also included the restrictions derived 
form the consistency condition (3.3) and the vacuum stability constraint (3.1), which are 
relevant for A^; > and A^; < 0, respectively. 

To discuss solutions for A^; shown in Fig. 7 it is useful to consider the following regions 
for {rrih.m^): 

• m^>mh 

In this region m<^ is large, so the t- and ^/-channel annihilation diagrams with (p 
exchange can be neglected. In this approximation (av) depends on A^, so there is no 
difference between left and right panels in Fig. 7. For A^; > 0, the constraint (6.11) 
is consistent with (3.3) in this mass range, however for negative A^; (3.1) requires 
< 1 TeV. 

For Aa; > the range of very heavy scalar masses ^1 — 8 TeV, corresponds to 
Xx ^ 1 — 10, which though large, remains below the unitarity and consistency limits. 
In this region of parameter space the singlets provide a substantial contribution to 
Smf^ that can ameliorate the hierarchy problem (see Fig. 4), and was utilized for this 
purpose in [16]. 



mh/2 <m^<mh 

In this mass range the interference between s- and t- and ^/-channel annihilation 
diagrams is relevant, so some differences between negative and positive A^^ solutions 
are visible. For large N and A^; > a small region is excluded by (3.3). 

^ mh/2 

In the vicinity of the resonance, the annihilation cross section is strongly enhanced, 
so that (6.11) can be satisfied only for small \Xx\. 
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Figure 7. The coupling A^; as a function of obtained from the requirement (6.11); for a the 
Higgs rrih = 130 GeV, N = 1^6 and 12 (red, blue and purple bands respectively), and A^^ > 
(left panel) or A^^ < (right panel). The blue areas in the left and right panels correspond to the 
regions disallowed by the consistency condition (3.3) and stability constraint (3.1) for A^^ = Stt, 
respectively. The thick black lines show the unitarity limit (3.2) saturated by \Xx\ = 47r. 



• m^< mh/2 

In this case the process (fip hh is kinematically forbidden, so only the 5-channel 
Higgs exchange graphs with no Higgs-bosons in the final state contribute; as a con- 
sequence {av) (X A^ and there is no difference between left and right panels Fig. 7 for 
this mass range. For A^^ > this mass region is excluded by (3.3), while for A^; < 
values of m^^ below m/j/2 are allowed depending on N. 

For scalar masses below ^ 1 GeV there are no CDM solutions. 
6.2 Feebly Interacting Dark Matter (FIDM) 

In order to thermalize with the SM particles, any other species of density n must have a 
thermalization rate F = n (av) larger than the expansion rate H. For relativistic singlets 
one can estimate F ^ A^r/(87r) while H ^ T'^/mpi, therefore equilibrium prevails for 
T < Teq where 

(6.12) 

If \x is very small (- 10"^) the species of DM being considered here does not equilibrium 
until very late times Teq ^ 1 GeV; for A^ < 7 x 10"^^, Teq < 2.7 and the singlets 
would not have equilibrated before the present epoch. 

In the following we will consider the case of feebly interacting DM (FIDM). We will 
assume that the DM number density / was negligible at the Big Bang: lim^^cx) f{T) = 
and solve the Boltzmann equation for = 1 (6.1) with this initial condition ^; this then 
determines the relic abundance: 

n^uh' = N^ = ^/ (6.13) 

Pcrit Pcrit 



^In practice, the Boltzmann equation (6.1) simplifies for this initial condition since the term on the 
right-hand side can be dropped: the DM particles do not annihilate when limT^oo f{T) = 0. 
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Figure 8. Top panels: FIDM solutions to the Boltzmann equation (6.1) with boundary condition 
/(T = oo) = 0, = 1 for two values of m^^ and rrih = 100 GeV (thin lines), 130 GeV (medium 
lines), 160 GeV (thick lines) and lA^^I = 10~^^ (bottom curves), 10~^^ (middle curves), 10~^ 
(top curves). The curve labeled /eq corresponds to the equilibrium distribution. Bottom panels: 
corresponding curves for (crv) /K in (6.1). 



where is the present photon temperature, and pcrU is the critical density. Following this 
procedure it should be remembered that the electroweak phase transition occurs at T^^w — 
300 GeV. Above this temperature, the only tree level contribution to the annihilation cross 
section comes from the XxH^Hif^ coupling, therefore 



a>3oo = |yi-^ (6.14) 

The case of FIDM was discussed first by McDonald in [31] and more recently also in 
[12] and [14]. In [31] the contribution from the process ^ H was included as a separate, 
and dominant, term in the Boltzmann equation. Here we obtain similar results following a 
different approach, where the presence of a non-zero width Th (see (6.5)) accounts for the 
creation and decay contributions of the Higgs boson. In other words, to avoid possibility 
of double counting we assume that the thermal generation of scalars is included in the 
processes DM + DM ^ SM + SM. 

The solutions of (6.1) for parameters relevant for insuring (6.11) are presented in top 
panels of Fig. 8. As seen from this figure, if A^; ^ 10~^ the scalars will be in equilibrium 
with the SM only for T < 100 MeV in agreement with the estimate (6.12); for smaller 
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Figure 9. Solutions to the Boltzmann equation for FIDM, for A/" = 1, 12 scalars and ruh = 130 GeV. 



Xx equilibrium is reached only for even lower temperatures. In order to understand the 
behavior of the solutions for / we plot in the bottom panels of Fig. 8 the factor {av) /K 
appearing in (6.1)) as a function of T. We see that {av) /K reaches its maximum for 
T ~ 20 — 30 GeV, and is strongly suppressed for T < 10 GeV. Because of this / becomes 
T-independent for T < 20 — 30 GeV, as observed in the top panels of the figure. 

In the case of non-equilibrium solutions of the Boltzmann equation, (6.11) is satisfied 
only if <C 1. For these small couplings we have {av) oc up to small corrections, in 
addition the term oc on the right-hand side of (6.1) is subdominant, so that to a good 
approximation / oc A^. It then follows that for FIDM the relic density depends on the 
DM parameters in the combination ^ A^m^^A^; this behavior is indeed observed in 

Fig. 9. In particular DM masses above 1 GeV requires A^; < 10~^^. 

If the scalar masses are ^ MeV and A^; > then (3.3) requires A^; <C 1, for example 
A^ < 1.65 X 10"'^ and A^: < 1.65 x 10"^ for = 100 MeV and = 1 GeV respectively. 
Note also that from Fig. 7 we conclude that there is no CDM solutions for A^; > and 
rrif^ < 70 GeV ; for smaller masses, A^^ is also very small and the situation reverts to 
the FIDM scenario. One should however remember that solutions with A^; < are also 
possible; in this case Fig. 7 implies there are no CDM solutions for < 2 GeV, then 
again the only solution is the FIDM shown in Fig. 9. 



7 Direct detection 

In this section we discuss constraints imposed on the model by searches for direct signals 
of DM particles scattering off nuclei. Even though the results of many experiments are 
available, here we will concentrate on the constraints obtained by the XENONIOO experi- 
ment [32] as they impose strongest limits on DM - nucleon scattering cross-section ctj^m—n 
in the mass range of our interest; at the end of this section we also comment on the recent 
results for the CRESST-II experiment [19]. 

The relevant scattering amplitude is described by the Feynman diagram in Fig. 10; the 
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corresponding cross section is 



C^DM-N 



4 9 



(7.1) 



where the sum runs over ah quark flavors rrin is the nucleon mass and are the nucleon 
form factors as deflned in [30]. 




Figure 10. The Feynman diagram for the elastic scattering of off a nucleon. 

In Fig. 11 we show the regions in the {m^^rrih) plane allowed by the XENONIOO 
limits [32] on ctdm-n- In the left panel the low region is magnifled to illustrate the 
two allowed bands. The one corresponding to the resonance region rrih ^ (see also 

Fig. 7); for these values the annihilation is amplifled to such an extent that (6.11) requires 
a very small coupling A^^ so that this region is allowed by the direct detection. Since the 
XENONIOO data is available for > 5 GeV therefore the vertical band of masses below 
5 GeV is also allowed. The right panel shows the large allowed area available for increasing 
^DM where the sensitivity of XENONIOO is reduced. 
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Figure 11. Constraint coming from the XenonlOO direct detection experiment combined with dark 
matter abundance in the case of CDM. Region forbidden for = 1 is marked blue, for = 6 blue 
and yellow, and N = 12 yellow, blue and red. The left panel magnifies the region of low scalar mass 
< 80 GeV, while the right one extends the mass range up to 600 GeV. 



Recently a positive result for a direct detection in elastic scattering of DM particles off 
nucleons was announced [19] by the CRESST-II collaboration; with two points selected by 
a maximum likelihood flt (Tab. 4 and Fig. 8 in [19]) Ml : {moM = 25.3 GeV, ctdm-n = 
11.6-10"^ pb} and M2 : {moM = 11-6 GeV, ctdm-n = 3.7-10"^ pb}. These cross sections 
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are far above lower limit from XENONIOO experiment [32], an issue that is yet to be 
resolved. As a complement to the implications derived above using the XENONIOO limits, 
we will also discuss the consequences the CRESST-II results would have in constraining 
our model. 

For each choice of ruh and N the requirement (6.11) determines a narrow allowed 
band in the (Xx^tti^) plane; an example is provided in Fig. 7 for the CDM case with 
rrih — 130 GeV and N — 1, 6, 12. It is natural to ask weather there exist ruh (in the region 
allowed by the present data) and N such that the corresponding band contains the most 
likely points Ml and M2 found by CRESST-II. The answer is contained in the top panels 
of Fig. 12, where the regions preferred by CRESST-II (within 2cr) are superimposed with 
predictions of our model for = 1, 2 and rrih = 110 — 140 GeV when satisfying (6.11) at 
the 3a level. As it is seen from the figure both regions are consistent with the model. For 
larger values of A^ our model is consistent only with the 2cr Ml contour for ^ 50 GeV. 
If confirmed, these data would provide a strong restriction on the scenario discussed in this 
paper; in particular a tightening of the allowed regions around the current central values 
would require A^ < 5. Moreover it has immediate and dramatic consequences for Higgs 
boson searches at Tevatron or LHC since the invisible h (p(p decay has a substantial 
branching ratio for the parameter regions consistent with Ml or M2. This is illustrated 
in the middle and bottom panels of Fig. 12, where the branching ratio is plotted against 
scalar mass m^. 

When < 0, the constraint (3.3) generates no restrictions and the branching ratio 
can be very large, BR}^^^^ ^ 1 (see middle panels of Figs. 12 and 13). In these figures 
we also plot maximum value of BR^^^^ calculated for |Aa^| that saturates the vacuum 
stability bound (3.1) on \Xx\ for a given quartic (f self coupling A^^ = Stt, 1,0.5 and 0.1, 
when nih = 130 GeV; the regions above corresponding dashed curves are excluded by this 
constraint. 

For Aa^ > the constraint (3.3) is important and excludes the region above the dotted 
curve shown in the bottom panels of Fig. 12 and 13. Within the allowed region the BR}^^^^ 
is smaller, but it can still reach ^0.8. 

Note, that if one anticipates the cutoff A to be below ^ 10^ GeV (as we do in this 
work) then the triviality upper limit for low A^ is roughly A^^ < 5, so the branching ratio 
could easily reach even 99% without any conflict with triviality. It is also worth mentioning 
here that if A^^ > 0, then small value of m^^ in the Ml and M2 regions together with the 
consistency condition (3.3) imply 0.06 < A^^ < 0.24, which is too small to ameliorate the SM 
fine tuning problem (see Sec. 5). On the other hand, if A^^ < 0, the vacuum stability bound 
(3.1) implies |Aa;| < 0.85 for \^ < Stt and 110 GeV < rrih < 141 GeV; then, as seen from 
Fig. 5, the allowed values for the cutoff A are smaller than for the SM. Therefore A^; < 
does not help solving the hierarchy problem regardless of the cosmological constraints. 

From the top panels of Figs. 12 and 13 one observes that for 1 < A^ < 12 low Higgs 
boson masses {rrih — HO — 120 GeV) are preferred by the Ml region of heavier (larger 
A^ favors smaller ruh). For the region of lighter (M2) roughly all Higgs boson masses 
fit data equally well, however only for A^ = 1 and 2. 

In conclusion, the CRESST-II positive result, if confirmed, would imply (within the 
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Figure 12. Top panels: singlet-nucleon cross section for = 1, 2 as a function of when (6.11) is 
obeyed, the color bands correspond = 110 GeV (green), 120 GeV (dark green), 130 GeV (red), 
140 GeV (dark red); also shown the regions Ml and M2 favored by the CRESST-II experiment 
[19] at 1 and 2 a level (black and blue closed curves, respectively), and the XENONIOO limit [32]. 
Middle panels: the branching ratio BR^^^^ for = 1,2 as a function of m^^ when A^^ < 0, for 
four values of rrih (color coded as above); the dashed lines show the value of BR^^^^ when A^^ 
saturates the stability bound (3.1) when the quartic (f self coupling equals A^^ = Stt, 1,0.5 and 0.1 
and for rrih = 130 GeV: regions consistent with both (6.11) and (3.1) lies on the colored bands and 
below the dashed curves. Lower panels: the branching ratio BR}^^^^ for A/" = 1, 2 as a function of 
m^^, when A^; > 0, and for four values of rrih (color coded as above); the dotted line show the value 
of BRh^00 when A^^ saturates the consistency bound (3.3): regions consistent with both (6.11) and 
(3.3) lie on the colored bands below the dotted line. 
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Figure 13. Similar as in Fig. 12, for N = 6 and N = 12. 



singlet extension of the SM) that the Higgs boson could decay mainly invisibly, escaping 
discovery at both Tevatron and LHC. The option of invisibly decaying Higgs boson have 
been considered previously in [33]. 

8 Self-interacting DM 

In spite of its many successes, the standard CDM cosmological model is also facing some 
difficulties when compared with recent observations. For instance, high-resolution N-body 
simulations have shown that the model generates a cusps in the DM density distribution 
in central regions of galaxies [34]. Another discrepancy concerns the number of subhalos 
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predicted by the model, which is at least factor of ten larger than the observed [35] number. 
Self-interacting DM was proposed by Spergel & Steinhardt [36] to cure these problems. 
Within this model, dark matter particles experience weak, non-dissipative collisions on 
scales of kpc to Mpc for typical galactic densities. This effect then generates a soft core in 
the inner regions of the dark halos, and it also ameliorates the overabundance of subhalos. 
The key requirement is that the mean free path of DM particles should be between 1 kpc 
and 1 Mpc in regions where the dark matter density is about 0.4 GeV/cm^. This model 
has attracted much attention [37]. 

In order to discuss this idea quantitatively we define the elastic scattering cross section 
per unit mass for DM consisting of N species. We then imagine a clump of DM particles 
scattering off another such clump: 



where z, j, etc. label the DM flavor (in our case, the components of (^), rrii denote the 
corresponding masses, rii the number density of z-th particles in the clump, and (Jij^kl 
the 2 on 2 cross sections. Because of the 0{N) symmetry the ni — n and rui — are 
flavor-independent, therefore 



(8.2) 



(^DM _ ^ ^f,k,l ^ij^ki _ an^ii + (TV - l)crii^22 + {N - l)ai2^12 
rriDM Nm^ 

We will define an auxiliary self-interaction cross section: 

^ ^ _ {mix, + 32myXl - 4m^(m^Ai + v'Xj + 2v'Xl)f 

where Ai is the coupling of the direct quartic interaction, A2 is the s-channel Higgs exchange 
coupling and A3 corresponds to the t- and u-channel Higgs exchange. Now the relevant 
cross sections are: 

c^ii^22 = cr^^i^Ax.O) (8.5) 
^12^12 = 2a^^(^,0,A,) (8.6) 

where the factor | in (8.6) and (8.5) comes from the combinatorics of the relevant diagrams 
and the factor of 2 in (8.6) corresponds to the absence of two identical particles in the final 
state. The relevant diagrams are shown in Fig. 14. 

In terms of the cross section per unit mass (8.2), the Spergel & Steinhardt hypothesis 
requires that 

2.05 • 10^ GeV-^ < < 2.57 • 10^ GeV"^ (8.7) 

rriDM 

In this section we will investigate whether this constraint is consistent with all the other 
restrictions that we have imposed on the multi-singlet extension of the SM. 
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Figure 14. Feynman diagrams contributing to the elastic (f(f scattering. 
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Figure 15. Regions in {X^p^m^p) space allowed by the Steinhard-Spergel hypothesis for = 1 
(blue/darker), = 12 (pink/brighter) region and the Higgs-boson mass fixed at rrih = 130 GeV. 
The left, middle, and right panels correspond to A^^ = 0, 0.1 and 1.0, respectively. The horizontal 
lines shows upper limits on A^^o allowed by the triviality condition for the cutoff A = 10 TeV and 
N = 1,12. 



For a given A^ and m/^, the condition (8.7) defines an allowed region in (A^;, A(^, TTiip) 
space; the projections onto the {m^^^X^^) plane for fixed \Xx\ — 0, 0.1, 1 with A^ = 1, 12 
and ruh — 130 GeV, are shown in Fig. 15. In the plots A^^ varies form up to Stt, the 
maximum value allowed by perturbative unitarity (3.2). It is clear that the condition (8.7) 
poses a very strong constraint on m^^ which is restricted to small values, generally well 
below 1 GeV for |Aa;| < 1 and A^ < 10 {m^ values ^ 100 GeV are allowed only for for 
|Aa;| ^ 10 and A^ > 100); similar results were obtained in other versions of scalar DM 
models [6, 31, 38]. 

If, in addition, one requires that the scale of physics beyond the singlet extension is 
above 10 TeV, then Landau poles are forbidden below this value for A ^. As a consequence, 
the limits on m^^ are much stronger, for example in the case of A^; = we obtain: m^^ < 
0.025 GeV (A^ = 1) or < 0.045 GeV (A^ = 12). 

As observed from Fig. 7 the values of rui^ that are required by the Spergel & Steinhardt 
condition (8.7) are not compatible with the CDM case. The only viable option is the 
FIDM. In this case, from Fig. 9 one can see that m^^ ^ 0.01 — 0.05 GeV corresponds to 

^In our derivation of aoM we neglected radiative corrections: because of the small allowed values of m^, 
we expect that the effects from the renormalization group evolution of the couplings will be small, and for 
this reason in this section we do not differentiate between Ac^ and its initial value (4.6). 
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Xx ^ 10~ — 10~ for = 1 — 12, and only the first panel in Fig. 15 is consistent with the 
DM abundance. Thus, within our model, the Spergel & Steinhardt hypothesis is consistent 
only with very light DM particles (m(^ < 0.01 — 0.05 GeV) that are very weakly coupled 
to the SM (A^ - 10-12 _ 10-10). 

9 Summary and conclusions 

We have considered an extension of the Standard Model that contains a set of N real scalar 
gauge singlets that transform as the fundamental representation of a global 0{N) sym- 
metry, under which all SM particles are singlets. This global symmetry remains unbroken, 
so that is a stable candidate for Dark Matter. The constraints on the model parameters 
implied by tree-level vacuum stability, unitarity and triviality were discussed. We have also 
investigated to what extent the presence of extra scalars could ameliorate the fine tuning 
of quadratic corrections to the Higgs boson mass. 

The restrictions imposed by abundance of DM were also presented; in particular we 
showed that the recent CRESST-II data for DM-nucleus scattering, if confirmed, would 
imply that for this model the standard Higgs boson decays predominantly into pairs of 
Dark Matter scalars. It that case discovery of the Higgs boson at LHC and Tevatron is 
impossible. For = 1 the most likely mass of the dark scalars lies in the range 15 GeV < 

< 50 GeV and BR{h 00) ^ 96%. If A^ > 2, the scalars have to be heavier: 
50 GeV < rUip < 60 GeV, and BR{h 00) ^ 98%. We have also shown that the Spergel- 
Steinhardt solution of the Dark Matter density cusp problem restricts the parameters to a 
Feebly Interacting Dark Matter region with - 10"^^ - 10"^° and scalar masses in the 
range - 0.01 - 0.05 GeV. 
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